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1. Introduction 

Brane interactions in the brane world have been proposed as the origin of inflation 
[H 121 E] 5 an epoch in the early universe that initiated the radiation-dominated big bang. 
In brane inflation |1] the inflaton is an open string mode, while the brane interaction 
comes from the exchange of closed string modes. A particularly simple scenario is the 
Dp-brane-anti-Dp-brane system. Although the potential between them seems too steep, 
it was proposed [Sj that the compactification effect in the special case of a hyper-cubic 
torus will flatten the inflaton potential for enough inflation. Let us call this the DD 
scenario. 

Recently, it was pointed out jHllZj that the Poisson equation for the inflaton potential 
$ in a compactifled manifold should include a background term (the so-called "jellium" 
term in solid state physics). As a result, the slow-roll parameter rj in the DD inflationary 
scenario (in the simplifled version where the stabilization of compactification moduli is 
independent of the inflaton) becomes 

V - -2/rf± (1) 

where d± is the number of dimensions perpendicular to the branes, d± = 9 — p < 6. 
Since we need at least Ne = 50 e-folds of inflation, and |?7| < l/N^, the DD scenario is 
not viable as an inflationary model. The impact of the jellium term is clearly important 
to the analysis of the inflationary properties in some of these scenarios. This was flrst 
studied in Ref. [S]. 

Here we would like to point out that the brane inflationary scenario where branes 
are at a small angle jHl E]] remains robust because the jellium term is much smaller 
in this scenario. In the {n,l), (n, — 1) wrapping scenario (which reduces to 2n parallel 
D4-branes after inflation), 

02 

ric:i (2) 

where 6 is the angle between the two branes (^ ~ 1/12 and n = 8 are reasonable values). 
When the jellium contribution to \rj\ is much less than 1/50, the slow-roll behavior of 
the inflaton is dictated by the other terms in the potential, in particular the quartic 
harmonic term as measured around the antipodal point, as studied in Ref. |10j. 

We would also like to point out that the DD scenario may still be possible under 
some special conditions. Whether that special condition on the background charge 
distribution is realized or not depends on the dynamics of moduli stabilization (the 
dilaton, the complex and Kahler structures of the compactifled manifold etc.), an issue 
that needs better understanding fn\ IT^ I7j. 

In section 2, we give a review of the Poisson equation in compact spaces. A simple 
ansatz of calculating the potential energy between two charges is gi$2 where $2 is the 
potential due to the the second charge q2, irrespective of whether the source charges in 
compact space add to zero or not. We shall justify this ansatz by showing that it is 
equivalent (up to a constant) to the potential energy between two charges (NS-NS or 
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RR) by directly integrating the potential energy density over the compact space. This 
equivalence is simply illustrated in the 1-dimensional case jHj- In higher dimensions, 
the Green's function in compact space may be represented in terms of Jacobi theta 
functions. In section 4, we generalize the method used in solid state physics to show 
how to obtain the Green's function numerically. This method should be useful in more 
realistic brane world models. For flat compact spaces, the Green's functions obtained 
by these two methods agree. In section 5, we apply the result to brane inflation in the 
cosmological context. As pointed out in Ref. [HllZj , the slow-roll parameter rj is too big for 
the brane-anti-brane scenario. On the other hand, the impact of the jellium term in the 
branes-at-small-angle scenario can be made negligibly small. In this sense, the branes- 
at-small-angle inflationary scenario remains robust. Section 6 contains discussion. 

2. Poisson Equation and Potential Energy in Compact Spaces 

In non-compact space, the determination of the potential energy (Coulombic or 
gravitational) between two charges is well known: we treat one of the charge qi to be a 
probe charge, and the potential energy is given by gi$2; where $2 is the potential due to 
charge q2 (even if |gi| > \q2\)- Here, we argue that this simple ansatz is equally applicable 
in compact spaces, where $2 includes the contribution of the jellium term. That is, the 
potential energy between two charged (Coulombic, Ramond-Ramond, gravitational or 
NS-NS) objects in a compact space will include the same quadratic component due to 
the jellium term, irrespective of whether the sum of the source charges in the compact 
space vanishes or not. 

Consider a compact manifold M {dM = 0). The Green's function is given by the 
Poisson equation. 



where 6 {r — r') is the cix-dimensional 6 function and V is the volume of the compactified 
manifold. Integrating both sides over the 0?^- dimensional compact space and using 
Stoke's theorem, we obtain Gauss's law. Since a compact space has no boundary, this 
volume integral over V^G vanishes, as does the integral over the RHS of Eq.Q. 

Consider the eigenfunctions ux (r) of the Laplacian operator on the compact 
manifold with eigenvalues A < 




(3) 




(4) 



The eigenfunctions satisfy the completeness relation: 




(5) 



A 



Integrating both sides of Eq.Q over M, we obtain {dv = d'^^r 





(6) 
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In a compact space there is a normalizable zero mode of the Laplacian, V^Mq (i") = 0, 
where Uq (r) = l/^A^. The Green's function can be written as: 

G(r,r') = V^^^^^i^ (7) 

Note that the zero mode is absent in G(r, r') in Eq.Q, thus avoiding an obvious 
divergence. Using Eq.® and we have 

/ dv G{r, r') = (8) 

J M 

Let $ be the static potential due to a given set of charges qi (at Vi), i = 1, 2, 

^r) = J2Q^G{r,r,) (9) 

i 

Using Eq.Q and Q, the total static energy due to this set of charges is 

V{rk) = ]- ! dviy^f = - [ dv ^$V2$ (10) 



Using Eq.® and Eq.(jni), we see that the last term in Eq- lfTUj) vanishes. So, irrespective 
of whether or not the net charge vanishes, 

\^(rfc) = — qiG{ri, Yj)qj + constant (11) 

i>j 

where the constant comes from the rj-independent self-interaction terms. When applied 
to the RR interaction, where qt = fii, we must have ^j/Xj = in our compact space. 
For the NS-NS interaction, qi are the brane tensions Tj, J^i'^i yNs-Nsi^^i) has 
the opposite sign since same sign charges attract. Putting them together we obtain, up 
to a constant 

V{rk) = + Y,{TiTj - /i./ij)G'(r„ r,) (12) 

i>j 

where the presence of the jellium term is encoded in G{vi, rj). In the two brane potential, 
the "Coulomb" {l/r'^~'^) term and the jellium (r^) term have the same effective strength, 
namely T1T2 — /ii/i2- This effective strength is very small for branes at a small angle. 
Since the constant term in the potential is linear in the brane tension T, V"/V can be 
made arbitrarily small. 

2.1. One- dimensional Example 

This well-known example (given in Ref. [H]) illustrates the above result. Consider a circle 
of circumference L. For < y < L, 

_ . . v"^ \v\ L . . 
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where the constant term may be determined by the ^-function regularization of the 
divergent Green's function G of the 1-dimensional lattice, 

oo 

2G{y) = \y\ + J2\y + ^M + \y-^L\ (14) 

n=l 

Notice that G{y) satisfies Eq.®. 

Let us put a charge gi = g at ?/i = and a second charge q2 = —q at y2 = y, where 
< y < L so that the total charge in the circle is zero. Using Eq. (fTT|) . one finds 

v(y) = «=Gfe) + # = ,'fM_l!^£W^^ (15) 



12 " V 2 2L 12J 12 
where the first term is the linear potential and the quadratic term is due to the jellium 
effect. The last term is the self-energy term and is independent of y. As we put the 
charges on top of each other, V^(0) = 0, as expected. 

We may also calculate the potential at x due to these two charges (+g and —q) : 

X 

(t)i{x) = q-{L - y) - qy + (t)Q 0<x<y (16) 
(p2{x) = — — + 00 y < X < L 

Note that (j){x) is continuous and piecewise linear in x, i.e., 0i(O) = 4>2{L) and 
4>i{y) = (p2{y)- Now we can integrate over the circle the energy density ((V0)^/2) 
to obtain the energy as a function of the relative position y of the charges. This exactly 
reproduces V{y) in Eq. (fT5|l . In short, (t){x,y) is the potential for a probe charge at a;, 
while V{y) is what we are interested in. 

For the case of two masses on a circle: nii at ?/i = and m2 at y2 = y, the potential 
at X (0 < X < L) due to these 2 masses is (f){x) = miG{x) + m2G{x — y), where there is 
a jellium contribution in each Green's function. From this, 

„ , , , mi + m2 Tn2y mi — mo , 
V0i(x) = L—^x+^+ 0<x<y (17) 

mi + m2 m2y nii + m2 ^ ^ ^ 

V(p2{x) = x + — y<x<L 

and, as expected, we have 

V{y) = -- [ dx{V<Pf = m,m2G{y) ^""^ + ""'^^ 



where the last constant term is the self-energy contribution. 



3. Green's Function in Arbitrary Dimension 

When two branes are moving in a compactified space, we need to calculate the inter- 
brane potential. This requires finding $(r) = — G(r) that satisfies Eq.(j2I). Our main 
goal is to determine the Green's function around the antipodal point (the point farthest 
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from the source). To simplify the problem, consider a square, flat d-dimensional torus 
of volume L'^. Using the identity 

oo 

e^'^'ds (A < 0) (19) 



L2A 

and the eigensolutions 



M„(x) = yZ^e^---/^ An = — (20) 



L2 



we can regularize 



p27rin-x/L 

GW= -Ea-Mx;;^;;^^ (21) 



poo d, 

/ ~ n e2™^"^/^-'"'"^)cis 



'2-c 



'2-d 



d 

e 



i-n»3(^. 



^■^1' 471-'^ 



ds 



Where 6*3 is a Jacobi theta function [12] 

^3 (^,e-M = ^6"''"'+'^" (22) 

^ ^ n 

The rapid convergence of the above expression makes the flat (i-torus especially simple 
to numerically simulate. 

When one considers compactification on a manifold (T^, K3, CY3) with a non-trivial 
metric, another numerical method suggests itself. Developed by Ewald and others 
this method is suited to spaces that may be represented as a periodic lattice (perhaps 
modulo some additional discrete symmetries). To illustrate the method, we proceed 
with a simple example, a square, fiat torus. Let us again write $ as in Eq.(|7)): 

/ N 1 expfik,- ■ r) , , 

^(r) = V E ,2 ' (23) 

where e ^Z'^ is a reciprocal lattice vector and kj its magnitude. This sum diverges 
for d > 2, since the number of terms grows as kj~^ (i.e., with the surface area of a 
(i-dimensional ball). Alternatively, one may treat the torus as an infinite lattice and 

sum over the source and its images at the lattice points: 

a 



where 



a = ^4^, d^2 (25) 

4:712 

As we saw in the one-dimensional case, this will lead to a divergence and the need to 
regularize it. The answer depends on how the lattice points are summed and how the 
regularization is carried out. 
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Ewald's method is to add and subtract a diffuse charge distribution around each 
lattice point. The added charge makes the monopole moment of each cell zero, aiding 
the convergence of the real lattice sum. This additional charge (minus a jellium term) is 
convergently subtracted away in reciprocal space, since a reciprocal space sum converges 
for non-singular charge distribution. The final result is independent of the diffuse charge 
distribution used to regulate the sums. Details on how to implement this method in two, 
three and four dimensions are included in the appendix. The numerical data generated 
by the Ewald method may be fit to an expression for the potential near the source of 
the form: 

(^) = ^ + ^ + "^i—clr) + constant, d > 2 (26) 

" ~ ^ (z) + iV + *i~.c(r) + constant, d = 2 
Notice that although individual terms in ^^'^\r) are not periodic, their sum is. The 
constant is fixed by requiring ^^^\r) to be independent of the added and subtracted 
charge. The harmonic piece satisfies the Laplace equation (when restricted to the cell 
containing the origin): 

V^$ilL = (27) 
For a hypercubic torus (with V = L'^) it has the general form 

= ^ + Afhf\r) + Afhfi.) + ... (28) 

V d 

where the hn^ are polynomials of order n with coefficients determined by Eq (j27j) . For 
example, for a hyper-cubic lattice with coordinates Xi measured from the source (podal 
point) and r = {xi,X2, ■■.,Xd) 



1 



- -m 

V d 



2 2 



(29) 



1 



hf\r) - -SI 

V d 



^2 2 2 , {d-2){d- 1) 6 _ id-2) ^ 4 2 



In two dimensions, terms of order 4m + 2 are not present in the harmonic piece. Beyond 
the sixth order terms in the hypercubic case, and for d > 2 for a rectangular lattice, 
there is more than one undetermined coefficient A^^] at each order. A two-dimensional 
rectangular lattice has only one parameter at any even order. At a given order n, 
the number of parameters for a hypercubic lattice reaches a maximum and becomes 
independent of dimension for d > ^. 

There is an expression similar to Eq. ()26|) for the potential near the antipodal 
point. This is the expression that is suitable for applications in brane inflation. With 
coordinates Zi now measured from the antipodal point (z = r — {L/2,L/2, ...)), $ has 
the form 

CWo...(-) = ^ + ^ E + B^'f^^'i-) + Blf'}^^-^ + ■■■ (30) 

i=l 
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Table 1. Constant, fourth order and sixth order coefRcients in potential near source. 



d± 


2 


3 


4 




-0.21 


-0.21 


-0.17 


A, 


0.12 


0.44 


0.34 


A, 


0.00 


0.0072 


3.05 



Table 2. Constant and coefficients of the fourth order and sixth order terms in the 
potential near the antipodal point . 





2 


3 


4 


Ca 


-0.055 


-0.064 


-0.070 




-0.62 


-0.53 


-2.20 




0.00 


0.0024 


-101.5 



The results for the coefficients in Eq. ()26j) and Eg. ()3()|1 were obtained by the method 
described below, and are summarized in Tables ^ and |21 The lattice spacing has been 
set to one. At least four terms in ^harm were kept in each dimension, but the accuracy 
of the numerical values could be improved by keeping more. In general, the convergence 
of Eq. is somewhat better than that of Eq. (jHOI) • 

It is easy to check that the integral of $ over a unit cell is zero, that is, $ satisfies 
Eq.®. In solid state physics, the Madelung constant is found by considering the 
potential due to both the positive and negative ions. Since the negative ions are found 
at the antipodal point in a simple cubic lattice, the Madelung constant is given by 
«m = Cs — Ca- In three dimensions, our value for the Madelung constant of a simple 
cubic lattice agrees with Ref . [T3] . The constants Ca in Table El agree with Eq. (j^H) 
evaluated at the antipodal point, and are four times those listed in Ref. 

For rectangular torus, there will be quadratic harmonic terms of the form zf — zj. 
Their impact on inflation is discussed in Ref. . The hypercubic way to sum the lattice 
generates only the harmonic terms. The numerical values of B^ is at least a factor of 
3 smaller than that given in Ref. |10j. This weakens the potential and improves the 
inflationary scenario. 

4. Application to Brane Inflationary Scenarios 

Let us consider a few brane inflationary scenarios, where the moduli stabilization effects 
are ignored. In order to flnd the potential of the inflaton as seen by a 4D observer we 
need to calculate the low-energy effective action for the brane system. The interaction 
between the branes due to the exchange of closed strings depends on their separation, so 
we will decompose the coordinates of the two (stacks of) branes into the center-of-mass 
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and the relative separation. 

Assuming that the branes wrap n and m times the volume V\\, the low-energy 
effective action is obtained by expanding the DBI action: 

~ (m + n)£|,rp + nrp / (F+^^J-d^y^d'^y^ + mr^ [ (F+%\d ^y^d^ y2 (31) 



where Tp is the brane tension 

r, = Mf+V(27r)^^?, (32) 

where Mg is the string scale and Qs the string coupling. The coordinates of the brane in 

the transverse directions are expressed as: 

m 

yi = ycM H : — yr (33) 

m + n 

Tl 

2/2 = ycM \ — yr (34) 

and substituting these into the expression for S'e// we obtain: 

Seff = '^P 2{mln) J j '^'^^^y^^'y^ = J d'^l^.^^^ij (35) 
The relationship between yr and ip is given by: 



V = yrJ , „a ^pM| (36) 



mn 
{m + n) 

where V\\ = J dP-^^ = 
4.1. The DD scenario 

The DD potential is given in Ref. 1101, with d± = 9 — p : 

Viy)=2r/l'-—^ (37) 



where 



K = SvrGio = ^,^o (3^ 



2Mf 

and (3 = 2a given in Eq. ()25|) . Measured relative to the antipodal point, the position 
of the anti-brane is given by z, which is simply y shifted by half the lattice size. The 
important pieces of the potential for the slow-roll condition on inflation are the constant 
term and the quadratic piece due to the jellium term : 

V{z) = 2Tpil-'-z'^^ (39) 

where V = V±. As pointed out in Ref.|ni [Z|, the relevant slow-roll parameter r] is given 
by: 
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where the derivative is taken with respect to the scalar field i/j that appears in the low- 
energy effective theory, Eq.ijSHI)- Since d_i_ < 6, the slow- roll condition is never satisfied 
in this case, that is, the branes will collide far too early for any significant infiation to 
take place. 

A priori, it is still possible that the stabilization dynamics of the extra dimensions 
has some unusual features that suppress rj and realize the condition required for the 
viability of the DD infiationary scenario. One such possibility utilizes a warped 
geometry to suppress the inter-brane attractive potential [Zj. 

Here, let us consider 

V^-l- = J2 - r.) - Fir - (L/2, L/2, ..)) (41) 

i 

n^) = ^^ + /(^) 

where /(z) is multi-periodic and consistency requires it to satisfy 

dvf{7) = (42) 



' M 

It is not hard to imagine that /(z) originates from the stabilization of the moduli in 
the extra dimensions. As one can easily see, DD infiationary scenario is viable if i^(z) 
vanishes at the antipodal point. To suppress the quadratic term in the infiaton potential, 
we need to decrease the value of -F(z) at the antipodal point so that 

|r/| ~ |2VF(0)/c/x| < l/Ne (43) 
Suppose, for a torus, measured with respect to the antipodal point, 

^ l-ncos(fc,2:,) ^ E(fc^%)' (44) 

so that F{Q) = 0. This implies that the infiaton potential does not have an anharmonic 
quadratic term around the antipodal point. Such a DD scenario will be able to give 
enough infiation to render the model viable. (In fact, one may choose F{z) to reduce 
the contribution of the quartic term to rj as well.) Among other factors, the form 
of F{z) depends on the dynamics of moduli stabilization, an issue that is not fully 
understood [m 1121 IZ|. Generically, we should consider -F(O) << 1/V as a fine-tuning. 
Since r] is generically around 1 for the DD system, we need a fine-tuning of 1 in 100 
on F{z) to suppress 1]. In more realistic constructions of string models, it will be very 
interesting to see how such a condition can be satisfied. 



4-2. Branes at a Small Angle 

Let us consider the simple scenario P where p = 4, dj_ = 4, but with different wrappings 
of the branes. The X'^,X^ torus, has sides £\\ and u£\\ and the branes wrap one-cycles 
on the torus, as shown in Figure 1. These branes are separated in the X^,X'^,X^,X^ 
directions by a distance y. 
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The Planck mass is given by 

= (45) 

We start with the example shown in Figure^where the wrapping numbers are (1, 1) and 
(1, —1). We shall consider small 6, so the angle between the branes is ^ ^ tan6' = 2u. 
The constant piece of the potential is given by: 

Vo = r/ii (2VI + - 2) ~ r/iiu^ (46) 

and the full potential is: 

V (y) = T^inu^ + harmonic (47) 

2d±V_i 2y^ 

where Q includes the contributions from both the NS-NS and RR sectors. 

Q = ( 1 - ^ - cos^^ ~ My/IQ-K ^ Myi2T: (48) 



27rsin^ V 2 

where the (1 — sin^ ^^/2) term comes from the NS-NS sector while the — cos 6 term comes 
from the RR sector. The sin 9 in the denominator comes from the length of the brane 
along the ^\\ direction, y measures the relative positions of the branes. With the above 
expressions for the potential and the Planck mass, we evaluate the potential l^(z) at 
the antipodal point. The relevant slow- roll parameter r] there becomes: 



For u ^ l/MsC-w — ctGUT — 1/25 (a reasonable choice, where acur is the standard model 
coupling at the GUT scale), the slow-roll condition is easily satisfied, and it is possible 
to obtain more than 60 e- foldings of inflation since ~ l/r/. 

A few comments are in order. As the jellium term contribution to 77 is very small, 
the number of e-foldings is dictated by the first non-zero harmonic term. Here it is 
the quartic term with strength given in Table El As we pointed out earlier, the 
obtained here are at least a factor of three smaller than those used in Ref . \Wi . Thus 
the inter-brane potential is weaker than Ref.jTHI uses, and so it is much easier to get 
sufficient inflation than they indicate. 
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The Madelung term Ca will shift the vacuum energy of the inflaton potential. As 
is clear from Eq. ()47|) and Table |21 this shift is positive, that is, it increases the vacuum 
energy term in the inflaton potential. Its contribution tends to decrease the magnitude 
of rj. 

After collision, the two branes at a small angle reduce to two parallel branes 
(horizontal in Fig.(P)) with zero vacuum energy. In an orientifold, the tension and 
RR charge of these two branes are canceled by the presence of orientifold planes. This 
implies that during inflation, the branes at angle will have a non-zero force with the 
remaining branes and orientifold planes. For example, the interaction between the 
(1, —1) D4-brane and a (1, 0) D4-brane is proportional to 

and its interaction with orientifold planes is also suppressed by the same factor of 16, but 
with opposite sign. These contribute a small correction to the interaction between the 
(1, —1) and the (1, 1) branes. One can always place the (1, l)-brane, the (1, — l)-brane, 
and the (1, 0)-branes at initial positions that are favorable to inflation. 

4-3. Other Branes- at- a- Small- Angle Inflationary Scenarios 

Next we consider branes wrapping the long dimension of the torus more than once. 
Suppose one branes has {n, 1) wrapping and the other has [n, —1) wrapping. After 
collision, we are left with 2n parallel branes. In this case the angle between the branes 
is ^ = 2u/n and the constant piece of the potential is: 

Vo = 2r/|| (Vn2 + - n) ^ Tpi^ — (51) 
The branes now intersect in points, so the charge Q is given by: 

Q^^^Ml^ (52) 

271 n-^ 271 n 
The relationship between and y becomes: 



iP = y^lTpi\\- (53) 



and the slow-roll parameter rj becomes: 

7] ~ (54) 

n 

For small u and large n, the slow-roll condition is easily satisfied, and it is possible to 
obtain many more than 60 e-foldings of inflation since ~ The other slow-roll 

parameter e = at the antipodal point. In the actual case, the slow-roll parameters 
are dictated by the quartic harmonic term where e is negligibly small. For a reasonable 
choice n = 8, we see that the number of e-foldings that can be obtained for inflation is 
further improved. In fact, the quadratic term is negligible in this case and the inflaton 
rolhng is dictated by the quartic harmonic term discussed in Ref . [TU] . 



Inter-brane Interactions 



13 



We may also consider branes with dimensionalities other than p = 4. In the p = Q, 
d± = 2 case the D6-branes span the X^, torus and wrap different cycles of the 
X^,X^ torus. They are localized in the X^,X^ torus and the interaction potential is 
logarithmic. The potential and the Planck mass are given by: 

Ml- " 



y [y) = 1 ~ ^ ^ ' 

where Q ~ Mf6^ /16tt. Again, in the small angle approximation, the relevant slow-roll 
parameter becomes: 

(M,r||)' 

The slow- roll parameter is again small in this case, and the end of the slow- roll is 
determined by the attractive logarithmic potential. The region yielding enough slow- 
roll is reasonably large and there is no need to fine-tune the initial conditions. 

It will be interesting to work out the situation of other brane inflationary scenarios 

dHUmilH]. 



rj ^ -^^^^ (56) 



5. Discussions 



In the simplified scenario discussed here, the DD inflationary scenario is not viable. On 
the other hand, the branes at a small angle scenario remains a viable model for inflation. 
Cosmic strings are generically produced towards the end of the brane inflationary epoch. 
Using the temperature fluctuation in the cosmic microwave background radiation to fix 
the superstring scale, the cosmic string tension arising from the brane recombination 
in the branes at a small angle inflation happens to be much bigger than that in the 
DD scenario ^01123. If branes-at-a-small-angle scenario is preferred, one consequence 
is that the cosmic string tension will be on the high side, up to values just below the 
present experimental bounds. This enhances the hope to test the brane inflationary 
scenario via the search of signatures of cosmic strings. 

As is clear from the analysis, the brane inflationary scenario depends on the 
dynamics of moduli stabilization. Presumably compactification moduli are stabilized 
by some string dynamics, or effective potential. The minimum of such an effective 
potential measures the cosmological constant. So understanding the moduli stabilization 
problem implies some understanding of the cosmological constant problem, or in a 
less ambitious framework, moduli stabilization must accommodate the smallness of the 
observed dark energy ^2] • Hopefully, brane inflation in string theory allows us to address 
this important issue. 

We thank Cliff Burgess, Nick Jones, Shamit Kachru, Juan Maldacena, Fernando 
Quevedo and Gary Shiu for many valuable discussions. This material is based upon 
work supported by the National Science Foundation under Grant No. PHY-0098631. 
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Appendix A. Ewald's Method 

Here we will detail Ewald's technique and extend it to suit our purpose. Consider the 
(i-dimensional torus as an infinite lattice and let $(r) satisfy 

V^$(r) = 5^-5(Arj) + l (A.l) 

j 

<l>(r) = $i(r) + $2(r) + constant 

V2$i(r) = J2 -Pi^^j) exp(-e2ArJ) + 1/V 

j 

j 

Where Arj = r — rj, and rj is the j*'* lattice vector. Here P(Arj) exp(— e^ArJ) can be 
thought of as a charge distribution inserted at the j^^ lattice site, normalized so that 

/ P{r) exp{-e^r^)dv = 1 (A.2) 

The charge distribution may extend outside the unit cell, but integrating one distribution 
over all space is identical to integrating all distributions over the unit cell. The effect 
of adding and subtracting a diffuse charge distribution is to regularize the summation 
over the lattice. The full potential $(r) is the sum of $2(1") (the real-space sum) and 
$i(r) (the reciprocal lattice sum). For a good choice of P(r) the sums are separately 
convergent, and together reproduce $(r) for the original point-charge distribution. 
Clearly, $(r) should be independent of P(r) and e, which are chosen to enhance the 
convergence of the sums. There is a range of choices of e (for our calculation 18 < e < 24) 
for which the results of Eq. ljA.lj) are independent of e. Other crucial conditions are listed 
and demonstrated in the next section. This procedure is well tested since it is widely 
used in solid state physics to compute the Madelung constant. 



Appendix A.l. Three dimensions 

In three dimensions, a simple choice is to add and subtract a Gaussian charge 
distribution at each lattice site. This gives P(r) = e^vr"^/^. Then we have 



<I,(3>(r) = 1 -erf(£Ary) 

An ^ 



-E 



2 477^ Ar 



where kj is the reciprocal lattice vector. Note that the following necessary conditions 
are satisfied: 

1) The terms in $1 converge faster than fcj^. 

2) $1 reduces to Eq(j2SI) when e ^ 00. 

3) The terms in $2 converge faster than Ar^^. 
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4) $2 goes to zero as e — oo. 

5) $2 reduces to J2j ^'^/^ ~^ ^■ 

Any other choice for P(r) must lead to potentials that satisfy these same limits for the 
method to work. 

In principle there is also an arbitrary constant of integration, but the requirement 
that $ be independent of e yields 

$(3)(r) = $f)(r) + $(=')(r)-^ (A.4) 

The e-dependent term comes from the missing zero mode in $i(r), as can be seen by 
differentiating with respect to e. The constant vanishes as e — oo, so that Eq. ()A.4|) 
recovers Eq. lf^ . This also satisfies Eq.®. In practice, a suitable choice of e increases 
the convergence of the calculation. In condensed matter a jellium is always used to 
make the lattice neutral, but a more general term can be easily incorporated into this 
technique by including the appropriate Fourier expansion in $i. 

The numerical potential generated by summing over the lattice can be fit 

to a potential of the form given in Eq. ()26|) 

= ^ + ^ + ^ + Af'hrix.) + Afhf{.^) + ... (A.5) 

where the quadratic piece comes from the jellium, and the constant is related to the 
Madelung energy of the lattice. The data generated by the Ewald method can also 
be used to fit the potential near the antipodal point. We expect this expansion to be 
harmonic, except for the term. With coordinates now measured from the center of 
each cube, this is 

^fL^M = ^ + ^ E ^? + Bf'hf'iz,) + B^h^'iz.) + . . . (A.6) 

i=l 

Of course, if the expression near the podal point is known exactly, Eq ()A.6|l can be 
obtained algebraically using zi = xi — ^, etc. Generally, though, the higher order 
coefficients Aio, A12, . . . will not be small, so it is best to use a numerical fit. 

Appendix A. 2. Four dimensions 

To implement the same procedure in dimensions other than three, we first need to choose 
a suitable form for P{r). We will need the Laplacian in d dimensions: 



vV(^) 



(A.7) 



The choice for P(r) is guided by the behavior of $1 and $2 in the limits mentioned 
below Eq ()A.3jl . In particular, each should converge at least faster than or 



since the number of points in the sums over the lattice grow as the {d — 1) power of r or 
k. Also, $2(1") should reduce to the usual Coloumb potential as r — > 0. Taking e cxo. 



Inter-brane Interactions 16 

which amounts to ehminating the added charge distributions, recovers Eq (j23p . In four 
dimensions, a simple choice is |13j : 



(A. 



The coefficients (a, b, ..) can be adjusted to give a single term in the charge distribution 
P{r) oc r^. The choice of how many higher powers of r to include in the above equation 
is a matter of taste and simply changes n. Consider a = |, 6 = 0. Then 

p6 2 

pW(r) = (A.9) 

„2 iA/i\ „,.v.^„-i, / ;.2 



-T«(4)/ A 1 >^ exp(-A;V4e2)exp(ikj -r) / fcj \ 
$(4)(r) = $S^)(r) + $5'^(r) 



8Ve2 

The constant in the last expression above guarantees that $ is independent of e. 
Appendix A. 3. Two dimensions 

In 2 dimensions, a useful form for the real-space part is 

$f (r) = —Yr exp(-eV'2) f 1 + aeV + 6eV=^ + ...^ rfr' (A.IO) 

To evaluate $2 in the Ewald sum, note that some convienent expressions for the 
exponential integral are: 

u ^-^ n ■ nl 

n=l 

= exp(-t) ^ , , ^ 1 (A.12) 

For the simplest choice a = 6 = 0, we have 

e2 



00 



P^^\r) = - (A.13) 

TT 

^(2)/ N 1 exp(-A;V4e2)exp(ikj- ■ r) 

^1 = V P 



$(2)(r) = $f)(r) + $f (r) 



4Ve2 

Extending the techniques from the previous three sections, we now have a well-defined 
way to evaluate the potential numerically for general dimension, charge configuration, 
and metric. 
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